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Abstract 

We have calculated the first and second order corrections to several deep inelastic sum 
rules which are due to heavy flavour contributions. A comparison is made with the existing 
perturbation series which has been computed up to third order for massless quarks only. 
In general it turns out that the effects of heavy quarks are very small except when Q ~ m 
or Q > m. Here Q and m denote the virtual mass of the vector boson and the mass of 
the heavy quark, respectively. For Q ^> m the radiative corrections reveal large logarithms 
of the type In Q 2 /m 2 which have to be absorbed in the running coupling constant so that 
the number of light flavours nj is enhanced by one unit. However this has to happen at 
much larger values of Q i.e. Q ~ 6.5 m than one usually assumes for the flavour thresholds 
which appear in the running coupling constant. An alternative description for the heavy 
flavour dependence of the running coupling constant in the context of the MOM-scheme is 
discussed. 
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The study of QCD sum rules, as represented by the first moments of the deep inelastic structure 
functions, has lead to a deeper insight of the behaviour of the perturbation series. This became 
possible after new techniques were invented to evaluate the Feynman integrals up to four-loop 
order. Examples of these techniques are infrared rearrangement |]J, integration by parts ||, 
and the .Reoperation ||. Also important was the appearance of new algebraic manipulation 
programs like FORM || which enables us to evaluate the complicated traces of the huge amount 
of Feynman graphs characteristic of higher order loop calculations. At this moment the sum 
rules computed up to third order in a s are represented by the first Bj0rken (polarized) sum 
rule |J, the second Bj0rken (unpolarized) sum rule |6|] and the Gross-Llewellyn Smith sum rule 
0. The perturbation series for these sum rules show a similar behaviour as is observed for 
other quantities which are calculated up to third order like e.g. the Z-boson and r-lepton decay 
widths (for a review of the literature see ||). Quantities computed up to a very high order in 
perturbation theory provide us with a very good tool to understand methods used in improved 
perturbation theory. Examples are the principle of minimal sensitivity (PMS ||) end the effective 
charge approach (ECH These methods were applied || to the above sum rules to obtain 

an estimate of the unknown order oq contribution. Another way to get the latter term is to 
use Pade-approximants as carried out in [|ll|] (for an estimate using renormalons see also fl2|). 
One of the remarkable results of these methods is that all estimates agree very well with each 
other. Apart from the theoretical interest there is also a practical one. Quanities which can be 
calculated up to a very high order in perturbation theory provide us with an excellent tool to 
measure the running coupling constant a s . Notice that in many cases the perturbation series is 
only known up to next-to-leading order (NLO) which means that, apart from some resummation 
of dominant terms, we have no control on the higher order corrections. An example of the 
determination of a s is given in [|13| where it is extracted via the polarized Bj0rken sum rule from 
the data obtained for the longitudinal structure function g±(x, Q 2 ). 

The order corrections to the sum rules mentioned above have been carried out in [|14| (the 
unpolarized Bj0rken sum rule) and ]I5| (the polarized Bj0rken sum rule and the Gross-Llewellyn 
Smith sum rule). In these calculations only massless quarks were considered but mass effects 
coming from the contribution of heavy quarks were omitted. The latter are important because 
apart from additional corrections the mass effects indicate when a heavy quark has to be treated 
as a massless or as a massive quark. This also indicates which number of light flavours nj 
has to be chosen in the perturbation series in particular for the running coupling constant at a 
given value of Q 2 . Here Q denotes the virtual mass of the intermediate vector boson in deep 
inelastic lepton hadron scattering. Before presenting the heavy flavour contributions we first 
give the definitions of the three aforementioned sum rules and the corresponding perturbation 
series corrected up to third order in a s . The polarized || and unpolarized Bj0rken || sum rules 
are defined by 
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Table 1: The coefficients in the MS-scheme of the perturbation series (|4]) corresponding to the three 
sum rules in Eqs. (^-(Bl). 
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Here # c denotes the Cabibbo angle and for the constant K(rif) we have quoted the values given 
by the flavour group SUpijif) for rif = 3,4, where nf represents the number of light flavours. 
The perturbation series of the above sum rules in the case of massless quarks can be written up 
to third order in a , as 
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where /3o and /3i stand for the first and second order contributions to the /3-function which are 
given by 
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The other coefficients ai, bi, Ci, which are computed in the MS-scheme in |TJ| and Ul5| , are given 
in table [I]. As has been mentioned above, the order contribution to Eq. is not known. 
However, there exist some estimates. Here we will adopt the results obtained from PMS given 
in p. They will be denoted by 
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where the coefficient a4 MS (3) is given in table [T|. It turns out that the other estimates originating 
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Figure 1: Forward Compton scattering graphs for heavy flavour production: W + d —>■ g + c . The 
down quark and the charm quark are indicated by a dashed line and a solid line, respectively. 



from ECH || and the Pade-technique |L1] are very close to the PMS value. Besides the sum 



rules above we also have the Adler sum rule [|16[ given by 

AF 2 (Q 2 ) = f 1 - [iff (ar, Q 2 ) - F?(x, Q 2 )} = K(n f ) 
Jo x 1 J 

K(3) = 2 + 2sin 2 fl c (SU F (3)) K(A) = 2 (SU F (4)) , (7) 

which holds in all orders of perturbation theory. Furthermore it does not receive higher twist 
contributions or mass corrections. The latter we have checked in our computations presented 
below. The coefficients in table [l| are only determined for massless quarks (see [14 , |15|). In the 



subsequent part of the paper we will show how the perturbation series is modified by including 
mass corrections due to heavy flavour contributions. 

In our calculations we assume that in addition to the gluon the proton only contains three light 
flavours given by the quarks u, d, s, including their anti-particles. The heavy quarks only show 
up in the final state. Since the sum rules presented above only involve non-singlet contributions 
the perturbation series for heavy flavour contributions in the case of neutral current interactions 
starts in order a 2 . However, for the charged current interaction we get already contributions on 
the Born level. Starting with the latter interaction AF\ in Eq. (0) and AF3 in Eq. ([|) are in 
lowest order given by the flavour excitation process 

d (d) + W -> c (c) . (8) 

In the process above we have only considered charm production because the other heavy quarks 
are heavily suppressed by the mixing angles occurring in the Kobayashi-Maskawa matrix. More- 
over the relevant values of Q 2 are so small that they are far below the thresholds of bottom and 
top production. Further notice that the integrals over the strange quark and the anti-strange 
quark densities cancel against each other in the computation of AFi and AF 3 so that the pro- 
cess in Eq. (ra) does not contribute to the sum rules when d is replaced by s. In NLO one gets 
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Figure 2: Forward Compton scattering graphs for heavy flavour production: V + q—>q' + H + H. 
The light quarks q and the heavy quarks H are indicated by a dashed line and a solid line, respectively. 



contributions from the virtual corrections to reaction (§) and the gluon bremsstrahlungs process 
(see Fig. ffl). The latter is given by 



d (d) + W -> c (c) + g . 



(9) 



The coefficient functions corresponding to the structure functions Fi (i = 1,2, 3, L) have been 
computed for processes @ and @ in [[L7| f\ for m d = and m c ^ (see also Notice 
that due to the non-vanishing mass of the charm quark one gets already a contribution to Fp 
on the Born level from reaction (|8]). If the quark in the initial state becomes massive the above 
processes also contribute to the neutral current interaction where W is repaced by Z or 7. The 
computation of the coefficient functions for the neutral current interaction where the masses of 



the initial and final state quarks are equal has been treated in |19| . Integration of the coefficient 
functions over the scaling variable x provides us with the result for the unpolarized Bj0rken sum 
rule 



.4 



Fi,(l) 



(Q 



2 2\ 

m 



Art 



1 



2 
? 



6 



77 Mn(l + 



sin 2 fl c , (10) 



and the Gross-Llewellyn Smith sum rule 



.4 



Fs,(l) 



(Q 2 ,m 2 



1 



3(1 + 



Air 



1 + e 



+ 2 



1+e 



sin d r 



respectively, where Cf denotes the colour factor which in QCD reads Cp = 4/3. Furthermore, 
we have also checked that the corrections to the Adler sum rule in Eq. (0) are zero as expected. 
The above expressions have to be added to the light quark contribution A V (3,Q 2 ) in Eq. (f|) 
to obtain the 0(a s ) mass corrections to the sum rules in Eqs. ©, (|) with K = K(4). In the 

2 Notice that A 2 in table I of [||| should be K A and not K A /2. 
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perturbation series above the following limits are of interest. When Q 2 is much larger than the 
mass of the charm quark, i.e. £ — > oo, then the corrections in Eqs. fliPP and flTT| ) tend to zero. 
This means that the sum rules are presented in a four light flavour scheme. However, when 
the mass of the charm quark becomes much larger than Q 2 , i.e. £ — > 0, then the heavy quark 
mass corrections are non-vanishing. After adding them to the massless quark result A r (3,Q 2 ) 
in Eq. (§) one can extract an overall factor K which turns out to be K(3) which is the value 
in the three light flavour scheme. This is expected because for infinite mass the heavy flavour 
disappears from the theory. Unfortunately this does not happen for the Adler sum rule in Eq. 
(0) because it is insensitive to the mass of the heavy flavours. The next process which shows 
up in neutral current as well as in charged current interactions is given by gluon splitting into a 
heavy quark anti-quark pair (see Fig. 0) 



q + V — q' + Q + Q with V = -y,Z,W. 



(12) 



The coefficient function for g 1; which is the same as for F 3 , has been calculated in ||20|| . Notice 
that the heavy quark loop contribution in Fig. to the gluon self-energy U(p 2 



m 



where p 

denotes the gluon momentum, has been renormalized in such a way that 11(0, m 2 ) = 0. This 
implies that heavy quarks are decoupled from the running coupling constant. The result for the 
polarized Bj0rken sum rule and the Gross-Llewellyn Smith sum rule becomes equal to 
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where Tf stands for the colour factor which in QCD is given by Tf = 1/2 (for Cp see below Eq. 
flllD). The coefficient function for F% can be derived from the ones obtained for the structure 
functions F2 and Fl which are presented in [2T|. They are obtained using the same renormal- 
ization condition for the heavy quark loop contribution to the gluon self-energy as given above. 
The result for the unpolarized Bj0rken sum rule is 
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Like for the flavour excitation mechanism (Eqs. (|8|), @) we have checked that the gluon splitting 
process does not contribute to the Adler sum rule (0). We are also interested in the asymptotic 
expansions of the expressions above. In the case the quark mass gets much larger than the 
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virtuality of the intermediate vector bosons we get 
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The expressions show that for infinite mass (£ — 0) the corresponding heavy flavour decouples 
from the radiative correction which is a consequence of the renormalization condition for the 
gluon self-energy in the graphs of Fig. |[ When the virtuality Q of the vector bosons is much 
larger than the mass of the heavy quark, which implies that the latter behaves like a light flavour, 
we obtain 
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The leading terms in the expressions above which are given by the constant and the logarithm 
ln£ can be predicted by the renormalization group. The general form up to order becomes 
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Table 2: The contributions to the sum rules originating from the light quarks A(nf), Eq. 
the charm excitation A^\ Eqs. ([TP]) , (|TlD , and gluon splitting into heavy quarks Ajj (H = c,b), 
Eqs. (0). For a comparison we also presented the asymptotic expressions A^ mp '^ 2 ' , Eq. 
in order a 2 . 



7 



where the coefficients aj,6j,Cj originate from the light quark contributions presented in table [I]. 
Notice that in the equation above we have already substituted the values of coefficients Pq, Pi, 
Eq. (|5|), appearing in the series expansion of the /3-function. The large logarithms of the type 
\nQ 2 /m 2 , which originate from all heavy quark loop insertions like in Fig. can be absorbed 
into the running coupling constant. This is equivalent to a redefinition given by 
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If we add the asymptotic expression (|i9|) for the heavy quark loop contributions to the perturba- 
tion series for the light quarks in Eq. (|) we obtain after substitution of a s (rif, /i 2 ) the following 
result 

A^ mp (Q 2 ,m 2 ) + A r (n f ,Q 2 ) = A r (n f + 1,Q 2 ) . (21) 

Therefore for Q 2 3> m 2 we obtain the expression of the perturbation series for massless flavours 
again but wherein now the number of light flavours is enhanced by one unit. The question is at 
which Q 2 this will happen. Here we will give the answer for the charm quark because the deep 
inelastic sum rules are studied in the region 2 < Q 2 < 100(GeV/c) 2 . 

In our analysis we will use the three-loop corrected running coupling constant which satisfies 



the matching conditions |22 



a. 



(n f , A n/ , /i 2 ) = a s (n f + l, A n/+1) /i 2 ) at /j = m nf . (22) 



If we choose A 3 = 397 MeV/c (MS-scheme) we get a s (3,/x§) = 0.375 for /j 2 = 2.5 (GeV/c) 2 . 
These values were obtained from a comparison of the polarized Bj0rken sum rule with the data 
carried out in Hl3| . Following the matching conditions in Eq. (p2|) one gets a s (5,M§) = 0.122 
(here A5 = 259 MeV/c) which lies a little bit above the world average of a s (5,M§) = 0.119. 
Nevertheless we will use a s (3, /z 2 ,) as a starting point. The values for the heavy flavour masses 
are chosen to be m c = 1.5 GeV/c 2 , m;, = 4.5 GeV/c 2 and m t = 173.8 GeV/c 2 . Further the 
number of light flavours in A r , Eq. (|j), and the running coupling constant is taken to be nj = 3 
irrespective of the value of Q 2 . In table || we have presented for Q 2 = 2.5,10,100 GeV/c 2 
the light and heavy flavour contributions to the perturbation series denoted by A(3) and Ajj 
(H = c,b), respectively. In the case of -4(3) , Eq. (|), we only consider the exact perturbation 
series corrected up to order and omitted the fourth order estimate 5A PMS (3), Eq. (^), which 
is listed separately in table |2|. The charm contribution, represented by the order a s corrected 
quantity A^\ is given by Eqs. ([l~0|), (|Tl~D which does not appear in the case of the polarized 
Bj0rken sum rule. The remaining heavy flavour contributions show up in order a 2 and they are 
represented in the table by A^ , 4^ (see Eqs. (0), (|T4|)). The top quark contribution is so 



(2) 

small that it is neglected. Besides the exact results for A H we have also made a comparison 
with the asymptotic expression given by the order a 2 contribution i^ ymp '' 2 ' (see Eq. ([19])) 
derived in the limit Q 2 ^> m 2 . From table ||| we infer that the heavy flavour contributions are 
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rather small even when compared with the estimate <L4 PMS (3). Only in the case of AFi(Q 2 ) 
at Q 2 = 2.5 GeV/c 2 the charm component is of the same size as the order a 2 estimate and it 
amounts to 0.017 which is about 2% of the light quark contribution given by A Fl (3) = 0.847. 
This effect can be wholly attributed to the charm excitation mechanism (see Eqs. (|8|), (Q)) 
represented by which also dominates AF 3 (Q 2 ). At larger values of Q 2 , decreases and it 
becomes of the same order of magnitude as A^f 1 which is due to the gluon splitting mechanism. 
Notice that the bottom quark contribution is always smaller than the charm quark component. 
The behaviour of the heavy quark contributions follows from their asymptotic behaviour at small 
and at large Q 2 , see e.g. Eqs. (|15|)-(|T8|). At increasing Q 2 the charm excitation contribution 
A^F 1 is decreasing whereas the gluon splitting part A^ becomes larger. At Q 2 = 100 GeV/c 2 , 
which is about Q = 6.5 m, the latter gets closer to its asymptotic expression 

^asym P ,(2)_ However, 

there is still a discrepancy between the exact and asymptotic expressions which in the case of 
the sum rules Agi(Q 2 ) and AF 3 (Q 2 ) amounts to 15 %. For AFi(Q 2 ) this is much worse and 
the difference between the exact and asymptotic expression is 28 % w.r.t. the exact one. In the 
case of the bottom quark one needs much larger values before A 
Q 2 > 1000 GeV/c 2 . In order to get A^ ymp,(2) 



A [ h ] within 1 



I (2) 

h 



^asymp,(2) occurs f or 



one needs the value Q > 25 m. 
Hence we can conclude that the large logarithmic terms given by ln(Q 2 /m 2 ) start to dominated 
the heavy flavour contribution for Q > 6.5 m which means that from this value onwards the 
heavy flavour behaves like a light quark. Therefore only for Q > 6.5 m the large logarithms have 
to be resummed as is explained below Eq. (|20|) which will lead to a rif + 1 flavour description. 

has to be changed into /i = 6.5 m nf . Using this 



m r . 



This means that the matching condition /i 
new matching condition we get, starting from a s (3, 2.5) = 0.375 as our experimental input value, 
the result a s (5,Mf) = 0.114. The latter is very close to the value obtained in fixed target deep 
inelastic scattering experiments given by ct s (5,M§) = 0.113 ||23|| . However, from the analysis 
above we think that the matching conditions as presented in Eq. (^2]) are rather artificial. There 
is no specific scale where nature suddenly jumps from an n^-flavour to an nj- + 1-flavour scheme. 
Moreover the relevant scale in deep inelastic scattering is q 2 = —Q 2 which is spacelike rather 
than p 2 > 4 ■ m 2 lf which is timelike. Here p denotes the gluon momentum in the graphs of Fig. Q 
Therefore, in principle all heavy flavour channels may contribute for spacelike processes which 
proceeds via the coefficient functions rather than through the running coupling constant. The 
decoupling of the heavy flavours from the perturbation series is then ruled by the Appelquist- 
Carazzone theorem pi}] . In order to get more continuity between the large and small Q 2 regions 
we substitute on the l.h.s. of Eq. (|2l"D at rif = 3 the coupling constant by 



„ t Q\ „MOM/, ,2\ 

a s (fj, ,3) =a s (/x ) 



„MOM 




+ Ui ) + 
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with 



Ui = T f £ 
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with m^ = m c , m^ = mj, m$ = m t 



(24) 



Here we are starting from a low input scale fiQ so that one is sure that at this value the per- 
turbation series is described by a three-flavour number scheme. The functions Hj are the order 
a\ contributions to the gluon self-energy which can be attributed to the heavy quark loop only. 
They are presented in |25| up to order a 2 . For /i 2 3> m 2 and m 2 lf 3> /z 2 , the expression in Eq. 
(|23| ) tends to its asymptotic result in Eq. (|20| ) provided one considers in the sum of Eq. ( p4[ ) one 
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flavour only. Finally one can resum the self-energies so that our new running coupling constant 
becomes 



MOM/ 2\ a s(H ;3) 



in 



The coupling constant above was proposed in the context of the momentum subtraction (MOM) 
scheme in [BfJ (for earlier work on this subject see also |2"Tf) where it also includes the resum- 



mation of the light quark contributions which we did not perform in this paper. The procedure 
outlined above guarantees that one gets a smooth transition between the low and large Q 2 re- 
gions. When e.g. /z 2 = Q 2 and Q 2 ^> m 2 the large logarithms ]n(Q 2 /m 2 l ) occurring in the 
functions U,i ( |24| ) cancel against the corresponding terms in the perturbations series of the heavy 
quark component AH(Q 2 ,m 2 lf ) of which the asymptotic expression is given by Eq. (|l~9"D. In 
this way one gets effectively a {rif + l)-flavour description. If Q 2 <C m 2 then Ui ~ and no 
large corrections appear in Ah(Q 2 , m 2 ) (decoupling of heavy quarks !) so that one gets the n/- 
flavour representation for the whole perturbation series. Following our approach and choosing 
//q = 2.5 GeV/c 2 in Eq. ( p4| ) we get a s (5, Mf ) = 0.117 which is closer to the LEP measurement. 
Finally we checked that the numbers in the table hardly change in passing from the rif = 3 
MS-scheme to the MOM-scheme. The maximal deviation is observed for Agi(Q 2 ) where at 
Q 2 = 100 GeV/c 2 the latter scheme leads to a decrease of A 9l (3) by 0.002 so that for this value 
of Q 2 the resummation effect is very small. Summarizing our findings we have calculated the 
heavy quark contribution to several deep inelastic sum rules. The corrections to the three-loop 
corrected perturbation series, computed for light quarks, are very small. Only at low Q 2 the 
correction in the case of the unpolarized Bj0rken sum rule AFi(Q 2 ) is noticeable where it is of 
the same order of magnitude as the order a^-estimate. The quark component of the sum rule 
attains its asymptotic value at much larger scales as given by the usual matching conditions. 
Matching at larger scales i.e. \i = 6.5 m nf leads to a smaller value of the running coupling 
constant at the Z-boson mass. The unnatural matching conditions which are characteristic of 
the MS-scheme can be replaced by expressing the perturbation series in the MOM-scheme. 
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